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1 .  INTRODUCTION 


The  prediction  of  radar  and  optical  sionatures  resulting  from 
atmospheric  nuclear  detonations  renui res  reliable  fireball  entrainment 
and  mixina  models.  Essential  to  these  predictions  is  a  sound  knowledqe  of 
the  temperature,  chemical  species,  and  velocity  fields  which  are  not  provided 
by  available  data  from  nuclear  tests.  The  objective  of  this  program  is 
to  investioatt.  turbulent  mixing  and  entrainment  of  air  in  fireballs  by 
means  of  a  subscale  laboratory  simulation  which  will  provide  a  data  base  to 
sunnort  the  development  and  evaluation  of  theoretical  fireball  models  and 
codes . 

The  development  of  firehalls  from  low  yield  ( < 1 00  kt)  bursts  at 
low  altitudes  (<30  km)  will  be  dominated  by  buovanc.v  forces.  The  initial 
stage  of  these  fireballs  resulting  from  large  enerqv  release  within  the 
atmosrhere  is  characterized  by  sot  erica!  symmetry,  the  domination  of 
radiative  enera.v  transport,  and  the  propagation  of  a  strong  spherical  shock 
wave.  On  the  order  of  a  second  after  detonation,  the  fireball  comes  to 
pressure  equilibrium  with  the  ambient  air  and  because  of  its  low  density, 
begins  to  rise.  As  the  fireball  rises,  circulation  is  aenerated  and  a 
toroid  is  formed.  The  series  of  laboratory  experiments  currently  being 
performed  are  desioned  to  simulate  these  "buoyant  rise"  fireballs  after 
pressure  enuilibration. 

This  renort  includes  descriptions  of  the  hardware  and  instrumentation 
added  to  the  test  facility  during  the  last  eight  months.  In  addition,  a 
new  ensemble  averaging  data  reduction  technique  is  described  in  detail  and 
reduced  interferometry,  shadowgraph,  and  particle  tracking  data  are  presented. 
Manv  aspects  of  the  current  investigation  were  discussed  in  detail  in 
Reference  1  and  are  not  repeated  in  this  report.  These  tonics  include  a 
support  analysis  and  a  full  description  of  the  high  pressure  test  facility 
and  holoaraohic  interferometer.  In  addition,  techniques  for  both  interferometry 
data  acquisition  and  Abel  inversion  data  reduction  were  described. 
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2.  EXPERIMENTAL  FACILITY 


Several  modifications  have  been  made  to  the  experimental  facility 
since  it  was  described  in  the  previous  Semi-Annual  Report  (Reference  1). 

The  more  notable  amona  these  involve  the  remote  controlled  burstinq  of 
the  helium- filled  soap  bubbles  and  the  addition  of  a  particle  seeding  and 
tracking  system.  A  nhotoqraph  of  the  nresent  overall  experimental  facility 
includino  the  hiah  pressure  test  tank,  holocamera,  and  other  support 
eouinment  and  diaonostic  instruments  appears  in  Figure  1. 

2.1  HELIUM  RELEASE  SYSTEM 

As  a  result  of  large  initial  asymmetries  present  in  vortices  generated 
bv  electrical  discharoe  bubble  breakinq,  a  mechanical  breaking  technique 
was  developed  which  reneatedlv  produces  a  symmetric  helium  release.  This 
technique  consists  of  remotely  dropnino  a  125  mil  stainless  steel  sphere 
from  an  electromagnet  mounted  inside  the  tank  throuqh  the  top  center  of 
the  bubble.  Tests  with  no  bubble  have  shown  that  no  more  than  1/16" 
dispersion  from  the  target  ever  exists  and  interferograms  taken  just 
milliseconds  after  arrival  of  the  sphere  at  all  test  pressures  indicate  that 
the  helium  release  is  indeed  symmetric.  A  schematic  of  the  new  bubble  bursting 
arranoement.  is  displayed  in  Figure  2. 

2.2  PARTICLE  TRACKING  SYSTEM 

Eouinment  providing  the  capability  of  obtaininq  multiple  exposure 
photographs  of  a  rising  vortex  with  its  initial  bubble  seeded  with  a 
particulate  tracer  has  recently  been  developed  and  incorporated  in  the 
test  facility.  It  consists  basically  of  two  parts:  an  additional  bubble 
formation  svstem,  and  a  camera/flash  system. 

2.2.1  Particulate-Seeded  Bubbles 

Several  different  types  of  flow  tracers  were  tested,  such  as  talcum 
powder,  fine  aluminum  flakes,  smoke,  and  hollow  enoxy  microspheres.  Smoke 
was  discarded  since  it  could  not  be  generated  at  the  high  operatina 
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pressures,  talcum  powder  could  not  be  photographed,  and  aluminum  flakes 
prohibited  reneatable  bubble  formation.  Consequently,  epoxv  microsnheres 
from  10  to  250  microns  in  diameter  were  incorporated  as  the  tracing 
"articulate. 

To  be  an  effective  fluid  tracer,  the  particles  used  must  possess 
relatively  low  inertias  and  hi  oh  drags  to  follow  the  streamlines  accuratel 
At  the  same  time,  they  must  be  sufficiently  large  to  become  visible  under 
intense  illumination.  Particle  trackino  is  particularly  valuable  to 
describe  flow  patterns  not  easily  accessible  to  other,  more  standard 
methods  of  measurement.  The  utility  of  this  method  depends  on  the 
accuracy  of  streamline  tracino  desired. 

Assuming  negligible  radial  velocities  and  net  external  forces,  and 
usino  Oseen's  approximation  for  CD(Re),  results  from  an  analysis  bv  Wright 
(Reference  2)  snow  that  the  path  dependence  of  a  particle  in  a  vortex  mav 
be  characterized  by  a  sinole  parameter, 


R0  / 6ra„D> 
■0  V  \  M '  / 


where  initial  position  of  particle  from  center  of  core  of  vortex 
VQ  =  initial  fluid  sreed 
a  =  radius  of  particle 
v  =  fluid  viscositv 

CnRe 

PQ  =  drag  parameter  =  =  1  +  0. 1875(Re)  -  0.0148(Re)2  +  ... 


(Reference  3) 


M'  =  apparent  mass  of  particle  =  i  7ra 3 (p  +  1  p) 

-J  P  L 


pp  =  density  of  particle 
p  =  fluid  density 


A  plot  of  this  particle  path  dependence  on  parameter  x  from 

Reference  2  is  presented  in  Finure  3.  It  is  evident  that  laroe  x  is 

-  •  LI  ~  '  o 

oesirabl e . 


For  the  Particles  used  in  this  study,  a  value  of  >  mav  easily  be 
computed.  An  example  is  "resented  to  illustrate  what  hannpns  to  a  100m 
diameter  narticle  within  a  vortex  havinq  a  core  radius  of  one  inch,  a 
circumferential  velocitv  of  one  ft/sec,  and  the  fluid  oredoninantl v  air 
at  one  atmosphere.  The  narticle  Reynolds  number  for  this  case  is 
calculated  to  be  2,  thus  the  draa  parameter  Do  usino  Oseen's  approximation 
is  1.33.  For  this  particular  case,  =  14  which  is  indicative  of  a  very 
oood  tracer.  Since  D0  increases  strongly  with  Re  whereas  M'  is  verv  weakl» 
dependent  on  Re,  the  rarticulate  becomes  a  better  tracer  in  the  hiqher 
pressure  reoimes. 

Two  features  of  the  bubble  blowing  arranaement  described  in 
Reference  1  prevented  it*  use  in  formina  bubbles  seeded  with  particulate. 
First,  the  hel  lum-rarticulate  orifice  in  the  bubble  tube  head  was  too  small 
and  rereatedl''  cloaned  with  soar  and  rarticulate.  Secondly,  with  a  laraer 
orifice,  one  could  not  obtain  a  soar  film  over  the  top  of  the  head  to 
enable  a  bubble  to  be  initiated.  Fo^  these  reasons,  a  new  bubble  tube  head 
and  mechanical  soar  amlier  we^e  develoned  and  incorporated.  A  schematic 
of  the  bubble  blowino  arranaement  with  the  rarticulate  mixer  appearing  in 
dashed  lines  is  also  Presented  as  part  of  Fiaure  2. 

2.2.2  Camera/Flash  System 

A  35mm  Nikon  F  with  motor  drive  havino  an  f5.6  macro  lens  was 
selected  as  the  ontimum  available  camera  to  acouire  the  narticle  trackinq 
ohotoararhs.  Kodak  Tri-X  Pan  (ASA  4DD )  was  chosen  for  film,  and  a  special 
processing  technique  usina  Acufine  developer  was  incorporated  which 
effectivolv  increased  the  sneed  of  the  film  to  ASA  1250,  noticeably 
enhancing  the  resolution  atove  that  obtainable  using  standard  developers. 

A  7 . 5kV  rower  supply  oreratinn  at  an  energy  level  of  about  50  joules 
was  developed  to  rapidly  charae  a  2.5uF  capacitor  bank  which  when  discharaed, 
iqnited  an  EG&G  FX98C-3  (3"  length,  70mm  diameter)  flash  lamp.  The 
capacitor  discharge  was  trigoered  by  a  standard  square  wave  signal  generator 
at  the  desired  frequency.  An  oscilloscope  sweep  triggered  by  the  release 
of  the  bubble  breaking  sphere  from  the  electromaonet  facilitated  accurate 
monitorina  of  the  flash,  enabling  the  time  at  which  each  exposure  was  taken 
after  bubble  breaking  to  be  accurately  known.  The  flash  was  masked  except 
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for  two  aligned  slits  to  provide  a  pulsed  vertical  sheet  of  light  directly 
above  the  bubble  tube  and  normal  to  the  line  of  sight  of  the  camera. 

This  illumination  system  permits  particulate  to  be  visible  in  only  a 
cross-sectional  slice  (less  than  1/2"  wide)  through  each  event  allowing 
detailed  visualization  of  the  internal  flow  structure. 
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Figure  1.  High  Pressure  Test  Facility,  Holocamera,  and  Support  Equipment 


Figure  2.  Schematic  Diagram  of  Bubble  Forming  and  Bursting  Systems 


3.  DATA  ACQUISITION  AND  INTERPRETATION 


In  this  section,  acouisition  and  interDretation  of  finite  frinqe 
interferograms  as  well  as  narticle  tracking  nhotogranhs  will  be  discussed. 

The  procedure  used  to  obtain  and  interpret  sets  of  infinite  fringe  inter- 
feroqrams  of  a  simulated  fireball  was  thoroughly  discussed  in  Reference  1 
and  will  not  be  repeated  here. 

3.1  ACOUISITION  AND  INTERPRETATION  OF  FINITE  FRINGE  INTERFEROGRAMS 

A  horizontal  view  finite  frinqe  interferoaram  mav  be  acouired  using 
the  same  procedure  as  an  infinite  fringe  interferogram  as  described  in 
Reference  1  with  one  additional  step.  This  step  consists  of  a  precise 
rotation  of  both  the  horizontal  beam  Mate  holder  and  its  adjacent  mirror 
between  acouisition  of  the  test  and  comparison  records  on  the  holociraphic 
plate. 

An  example  of  a  finite  frinoe  i nterferooram  is  presented  in  Figure  4. 

In  this  tvpe  of  interferogram,  the  test  disturbance  causes  a  displacement 
or  shift  in  the  frinoes  relative  to  the  undisturbed  portion  of  the  test 
section.  The  magnitude  of  the  fringe  shift  is  directly  relatable  to  fringe 
number  in  infinite  fringe  interferograms  and  therefore  to  an  integrated 
change  in  density  or  concentration  along  the  optical  oath. 

The  bulk  of  the  interferometry  data  to  date  have  been  obtained  from 
infinite  frinoe  interferograms  at  conditions  such  that  a  near  optimum 
number  of  frinaes  were  present  from  a  data  interpretation  and  reduction 
standpoint.  Generally  as  the  tank  pressure  was  increased,  the  vortices  were 
studied  at  a  later  time  in  their  development,  hence  at  a  higher  rise  position. 
Present  emohasis  on  interferometry  data  acquisition  (see  Section  6) 
focuses  on  conditions  at  which  either  too  many  or  too  few  fringes  would 
exist  for  interpretation  and  data  analysis  if  the  infinite  frinoe 
technique  v/ere  used. 

Even  using  the  finite  fringe  arranaement,  there  still  exist  frinqe 
spacing  limitations.  To  vield  nonambipuous  data,  the  fringe  shift  in  a 
given  lenath  of  an  event  cannot  exceed  the  number  of  fringes  laid  on  the 
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undisturbed  Portion  of  the  interferogram  over  that  same  length,  A  good 
pxamnle  which  illustrates  a  case  where  the  fringe  shift  per  unit  lenath  is 
area  ter  than  the  fringe  spacing  per  unit  length  is  an  i nterferoaram  of  an 
unbroken  bubble  where  very  large  frinae  gradients  are  present  near  the 
edae  of  the  bubble  (Fioure  5). 

Experiments  have  recently  been  completed  which  define  the  operating 
rarae  for  the  existing  finite  fringe  interferometry  apparatus.  The 
maximum  fringe  spacing  is  about  2  fringes/inch  which  will  be  used 
for  low  pressure,  high  rise  tests,  while  the  minimum  fringe  snacino  is 
about  20  frinoes/inch,  commensurate  with  high  pressure,  low  rise  event-. 

3.2  ACQUISITION  AND  INTERPRETATION  OF  PARTICLE  TRACKING  DATA 

The  procedure  used  to  acouire  nhotnnraphs  of  particulate  seeded 
vortices  is  verv  straightforward.  After  blowing  the  bubble  and  waiting 
a  few  seconds  to  allow  any  internal  motion  to  damn  out,  the  metal  sphere  is 
dropped  from  the  electromagnet  simultaneous  with  the  initiation  of  the  flesh 
1 iahti no  s vs  tern.  An  oscilloscope  sweep  is  triggered  by  the  electromagnet 
release  and  records  the  flashes  as  a  series  of  spikes  usino  the  output  of 
a  nhotodinde  mounted  inside  the  tank.  After  the  vortex  is  out  of  view,  the 
flash  is  shut  off  and  the  film  advanced  for  the  next  event. 

Interpretation  of  particle  tracking  data  in  this  sense  is  limited  to 
detailed  internal  structure,  and  trajectory  and  growth  observations  of  the 
rising  fireballs.  These  Particulate,  shown  in  Section  2.2.1  to  accurately 
follow  the  motion  of  the  fluid,  adenuately  illustrate  whether  or  not  a  vortex 
possesses  a  laminar,  wran-un  "jelly-roll"  structure,  or  a  turbulent,  more 
thoroughly  mixed  structure.  Unfortunately,  since  the  vortex  appears  many 
times  in  a  single  multiple  exposure  nhotoaraph,  the  constraint  of  not 
havino  vortices  overlap  together  with  the  difficulty  of  identifying 
individual  [articles  Precludes  a  direct  measurement  of  the  velocity  field. 
Acouisition  of  velocity  data  (both  mean  and  turbulent)  at  a  space  fixed 
point  with  a  laser  Doppler  velocimeter  (LDV)  is  discussed  in  Section  6. 
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Figure  4.  Typical  Eight  Atmosphere  Finite  Fringe  Interferogram 


Figure  5. 


Finite  Fringe  Interferogram  of  Unbroken  Bubble 
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4.  DATA  REDUCTION 


The  purpose  of  this  section  is  to  describe  in  detail  techniques  for 
reducing  holographic  interferometry  and  shadowgraph  data  as  they  apply  to 
the  present  study.  ^  y 


4.1  H0L03RAPHIC  INFINITE  FRINGE  INTERFEROMETRY  DATA 

A  data  reduction  technioue  has  been  developed  (Reference  1),  which 
wil  allow  mean  density  and  concentration  (i.e.,  mass  fraction)  profiles 
o  e  calculated  for  a  bimolecular  mixing  flow  given  a  fringe  distribution 
as  a  function  of  radius  for  an  axisymmetric  event.  An  examination  of 
ypical  infinite  fringe  interferometry  data  (Figures  9,  10.  14,  and  15) 
reveals  that  individual  events  are  in  general,  asymmetric.  It  becomes 
obvious  that  an  Abel  inversion  of  the  fringe  shift  eouation 


S  '  f  f  (n  - 


where  S  =  fringe  number  or  shift 

A  =  wavelength  of  light  in  vacuum 
L  =  integration  path  length 
n  =  index  of  refraction 

ds  =  differential  path  length  along  the  light  ray 
subscript 

observetMn  the^es^re6""1611  1 "  c°mPf son  scene  and  also 
the  test  scene  as  a  reference  condition 

is  not  valid  when  applied  to  individual  interferograms  such  as  these 
To  obtain  an  axisymmetric  event,  some  type  of  average  of  these  fringe 
istnbutions  must  be  obtained.  One  could  take  an  average  in  such  a  way 
to  yield  a  probability-type  result  for  density  at  a  specified  location  in 
space,  i.e.,  a  tank-fixed  approach.  This  would  result  in  a  very  smeared 
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out*  vortex  wM'ch  in  no  wav  would  resemble  a  typical  vortex.  Instead,  an 
amroach  was  developed  whereby  a  coordinate  system  would  be  established  for 
pach  individual  vortex  based  on  the  location  of  three  points  common  to 
all  vortices,  but  located  in  general  differently  from  event  to  event. 

This  approach  retains  the  important  physical  features  of  the  flow  and 
produces  a  symmetric  vortex  which  is  a  representative  average  vortex 
rather  than  a  smeared  out  vortex  average. 


4.1.1  Establishment  of  Coordinates  on  Individual  Interferograms 

An  interferoaram  is  vitally  examined  to  determine  the  position  of 
the  saddle  point  bv  counting  fringes  inward  from  the  edge  of  the  vortex. 

In  this  connotation,  the  saddle  point  is  defined  as  the  "oint  where  the 
frinoe  number  (see  Section  4.1.2)  simultaneously  attains  a  maximum  value  in 
an  axial  direction  and  a  minimum  value  in  the  radial  direction.  Two 
fringe  maxima  positions  (generally  having  different  fringe  numbers)  are 
determined  in  the  same  manner  and  an  x-axis  is  drawn  throuah  them.  The 
y-axis  is  then  drawn  orthogonal  to  the  x-axis  and  throuoh  the  saddle 
point.  The  print  where  the  x  and  y-axes  intersect  becomes  the  vortex 
coordinate  orinin,  or  for  the  sake  of  brevity,  simnlv  the  origin. 

Reference  to  Figure  6  will  show  that: 


x^  =  distance  from  origin  to  first  dark  frinqe  on  left  side 

x^  =  distance  from  origin  to  fringe  maximum  on  left  side 

x^  =  distance  from  origin  to  frinoe  maximum  on  right  side 

Xp£  =  distance  from  origin  to  first  dark  fringe  on  rioht  side 

v  =  distance  from  origin  to  saddle  point 

Tank-fixed  parameters  of  rurely  statistical  interest  are: 


n  =  rotation  of  x-axis  from  true  horizontal 

=  vertical  rise  distance  of  origin  from  ton  of  bubble  tube 
o  J 

x  =  horizontal  distance  from  origin  to  centpr  of  bubble  tube 
o  3 


*In  the  sense  that  the  averaged  event  would  be  much  larger  than  any  individual 
vortex  and  the  important  edge  gradients  would  be  reduced. 


The  center  of  each  fringe  is  located  and  numbered  where  it  intersects 
the  x-axis  defined  in  the  previous  section.  The  numbering  starts  with 
S  =  1/2  for  the  outermost  dark  fringe  (corresponding  to  the  half-wavelength 
shift  in  ontical  nath  discussed  in  Section  4.2  of  Reference  1),  1  for 
the  next  inward  lioht  fringe,  and  so  on,  increasing  by  1/2  each  time  a  new 
fringe  is  encountered.  Care  must  be  taken  to  give  the  same  fringe  the 
same  number  when  it  intersects  the  x-axis  more  than  once.  What  results  is 
m  ordered  pairs  (x...  S..)  j  =  1,  2,  ...  m  for  every  ith  interferogram, 
i  =  •••  n-  The  next  section  describes  how  an  average  fringe 

distribution  is  obtained  over  n  interferograms . 

4.1.3  Data  Averaoin in 

In  the  nrocess  of  establishing  an  average  vortex,  it  is  essential  to 
nreserve  the  important  physical  features  by  normalizing  the  characteristic 
dimensions  of  each  individual  vortex.  The  dimensions  selected  for  this 
purpose  are  the  locations  of  the  fringe  maxima  and  the  left  and  riaht 
edges  relative  to  the  origin. 

To  transform  each  event  into  a  system  of  coordinates  whereby 
maxima  collanse  onto  the  mean  maxima  and  edges  collaDse  onto  the  mean  edges, 
the  General  linear  transformation  x .  .  =  ax^.  +  b  is  used.  First,  the 
averane  locations  of  the  edges  and  the  maxima  are  computed  using 

XLE  =  n  £  XIE. 
i  =  l  1 

XLM  =  n  xiy. 
i  =  l  1 

XRM  =  n  5^  XRM . 
i  =  l  1 

-  1  n 

XRE  =  n  H  XRE . 
i  =  l  i 

The  event  is  then  divided  into  four  regions  as  shown  on  the  next  page  and 
the  normalization  is  accomplished  separately  within  each  region. 
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I  I 

I  I 

j Region  1  j  Region  2 


x 


LM  • 


Region  3 


0 


I 

I 

I 

xRMj 


Region  4 


x 


The  constants  for  a  and  b  in  the  oeneral  linear  transformation  are 
determined  by  matching  the  endnoints  within  each  region  to  the  averaoed 
enauoints.  For  example,  in  region  l ,  at 


\i  =  XLF .  ’  xi,i  =  xLEi;  and’at  xi,i  =  xlm.’  xii 


=  x 


LM. 


For  the  pth  data  n0int  on  the  itn  interferooram,  the  new  coordinates 
In  region  1 , 


.  th 


Xii  arc: 


XLE .  --  xii  1  XLM. 
i  -\ 


x .  .  = 


XLE  '  XLM 


l.?  x 


Lh  LNi 


x...  xij  +  XLM 


In  reoion  2, 


XLE  "  XLM 
XLE.  "  XLM.  XLMi 


XLM.  1  xi  i  1  0 


r. .  =  LM 

x|_m,  i  J 


In  reaion  3, 


°-Xi,i  ixRHj 


7  .  RM 

11  XRM -  XlJ' 

l 
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In  region  4, 


< 


< 


XRE  '  XRM  + 

XRE .  “  XRM.  XlVi  + 
i  i 


Once  alMx.j,  S..)  have  been  transformed  to  (x..,  S..),  there  exist  n 
Sii  vs'  xi^  nroflles  with  all  edges  occurring  at  and  x^,  all 
maxima  at  x^  and  xRH>  and  the  central  minimum  at  x  =  0.  These  profiles 
mav  be  averaoed  with  respect  to  S  at  constant  x  to  yield  a  mean  S  vs.  x 
profile.  It  becomes  readilv  anrarent  that  every  ith  distribution  will 


not  have  discrete  S  at  the  same  x.  What  is  done  is  to  choose  a  convenient 

xk  t0  be  used  for  each  i  event  (sa.v  x”k  =  ~  xLE,  k  =  1 ,  2,  . . .  L)  and 

1  i nearly  interpolate  between  (x...  S..)  and  (x..+],  S.,+1). 


The  innut  fringe  numbers  to  the  Abel  inversion  daiff  reduction  routine  for 

an  axisymmetric  averaqe  vortex  are  simplv  Si ,  =  -  T  S  for  k  -  l  ?  . 

_  -  k  n  ^Jik  or  K  '  '*  •••L- 

The  S  vs.  x  for  both  the  left  and  right  sides  are  computed,  and  provided 

that  the  number  of  events  taken  is  sufficiently  large,  should  for  all 

Practical  purposes  coincide.  S  is  then  extrapolated  smoothly  from  1/2  to 

zero  to  provide  an  outer  radius  of  the  vortex,  required  in  the  data  reduction 
procedure. 


The  number  of  events  n  required  to  form  an  acceptable  averaoed 
event  reouires  a  compromise  between  the  degree  of  asymmetry  in  the  averaged 
event  and  effort  reouired  to  acouire  and  reduce  the  data.  One  way  to 
investigate  how  many  events  are  required  is  to  Plot  parameters  indicative 
of  asymmetry  vs^  number  of  events.  Jhe  method  used  in  this  study  is  to 
calculate  x^,  xRE,  x^,  xR^,  S^,  SR^  for  different  numbers  of  events  used 
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andjlot  the  differences  in  sides  as  a  function  of  number  of  events- 
(^)  (^\ 

x  /  /  \c"  / vs'  n‘  ^en  these  parameters  converge 

E  M  jm 


E  AM 

sufficiently  close  to  zero,  a  satisfactory  nunber  of  events  have  been 
included.  Analysis  of  data  obtained  at  both  laminar  and  turbulent 
conditions  indicates  that  a  data  set  of  20  to  25  events  is  generally 
sufficient  to  satisfy  this  symmetry  requirement. 

4 • 1  ■ 4  Ejn^ebajJJjaj(e_ Interferometry  Data  Reduction 

In  addition  to  mean  radial  helium  concentrations  through  the  core 
of  the  vortex,  radial  helium  concentrations  at  various  oositions  through 
the  fireball  wake  may  be  extracted  from  the  interferometry  data.  One  is 
confronted  again  as  to  what  tyne  of  an  averaging  approach  should  be  used 
ano  commensurate  with  the  previous  sections,  a  free  coordinate  system  and 
normalization  with  respect  to  the  wake  edges  were  selected. 

The  coordinates  are  established  in  the  following  manner  (refer  to 
Figure  7):  An  axis  approximating  the  center  of  the  wake  is  determined 
by  eye  and  drawn  longitudinally  through  the  core  of  the  vortex.  The 
wake  coordinate  origin  becomes  defined  by  a  line  drawn  normal  to  the 
wake  axis  through  the  vortex  coordinate  origin  (where  it  was  defined  in 
Section  4.1.1).  All  cuts  through  the  wake  are  parallel  to  this  line  (normal 
to  the  wake  axis)  and  their  locations  are  based  on  fractions  of  vortex 
diameters  below  the  wake  origin. 

Fringe  numbers  along  these  cuts  are  recorded  in  the  same  manner  as 
was  described  in  Section  4.1.2.  The  fringe  profiles  for  the  same  z  /D 
cuts  are  averaged  over  n  interferograms  using  the  edges  and  center  as 
normalizing  parameters,  and  when  n  is  sufficiently  large,  a  svmmetric 
averaged  fringe  profiles  results.  Since  the  axial  synrotry  condition 
is  once  again  satisfied,  the  previously  discussed  Abel  inversion  of  the 
fringe  shift  equation  may  be  applied  to  calculate  radial  helium  concentrations. 
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4.2  HOLOGRAPHIC  FINITE  FRINGE  INTERFEROMETRY  DATA 

The  Abel  inversion  technique  mentioned  in  Section  4.1  applies 
eoually  well  to  finite  fringe  inter. Programs,  but  the  input  S(x)  are 
obtained  in  a  different  manner. 

The  procedure  used  to  reduce  finite  frinoe  interferoqrams  for  a 
horizontal  (in  same  direction  as  parallel  fringes  in  undisturbed  flow 
reaion)  cut  is  as  follows:  First,  by  inspection,  the  fringe  which  has  the 
maximum  shift  is  determined.  A  tracing  is  then  made  followino  this 
frione  across  the  disturbed  region.  The  fringe  snacina  can  be  measured  so 
the  vertical  displacement  of  this  fringe  from  the  undisturbed  may  be 
expressed  in  number  of  fringes  of  shift.  This  continuous  fringe  shift  as 
a  function  of  x  is  identical  to  the  discrete  fringe  number  as  a  function 
of  x  for  infinite  fringe  interferoarams . 

As  lonq  as  each  fringe  may  be  conti nuouslv  followed,  it  is  very 
straightforward  to  extract  several  fringe  distributions  through  each  event. 
However,  it  shruld  be  noted  that  this  simple  data  reduction  is  based  on 
horizontal  cuts  throunh  events  which  may  not  be  aligned  horizontally  (see 
Section  4.1).  Conseouently,  the  flexibility  inherent  in  the  present  infinite 
fringe  data  reduction  procedure  such  as  inclining  the  coordinate  axes 
through  the  maxima  is  not  possible  with  this  simple  finite  fringe  inter¬ 
pretation. 

It  is  Possible  to  reduce  data  along  a  line  not  parallel  to  the 
undisturbed  horizontal  fringes.  Through  a  very  thorough  inspection,  the 
two  frinoe  maxima  may  be  located  (if  the  event  is  skewed  these  maxima  will 
in  general  correspond  to  different  frinqes)  and  an  x-axis  can  be  drawn. 

Fringe  shifts  could  then  be  measured  on  a  fringe  bv  fringe  basis,  and  the 
data  reduction  would  proceed  as  previously  outlined  for  infinite  fringe 
data  in  Section  4.1.1. 

4.3  SHADOWGRAPH  DATA 

Shadowgraph  movies  are  utilized  in  this  investigation  to  study  gross 
motions  of  the  risinn  vortex  such  as  size  and  rise.  This  information  cannot 
be  obtained  with  interferometry  since  that  diagnostic  is  limited  to  obtaining 
data  at  one  instant  of  time  in  each  event.  After  acouisition  and 
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nrocessinq,  each  shadowgraph  film  is  projected  on  a  larae  sheet  of 
naner  where  the  outline  of  the  vortex  may  be  traced  frame  by  frame  until 
the  vortex  disappears  from  the  field  of  view.  Since  only  gross  features 
arnear  on  a  shadowgraph  frame  (chance  of  illumination  proportional  to  the 
second  derivative  of  density  rather  than  directly  with  density  as  is  the 
case  with  interferometry,  ?UCh  distinct  features  as  saddle  ooint  and 
maxima  cannot  be  distinguished,  but  the  outline  of  the  edae  can  be 
accurately  identified.  While  viewina  the  shadowgraph  movies,  one  maw 
nuite  accurate!*  visually  establish  a  radial  axis  through  the  center  of 
each  event.  The  data  averaoinq  procedure  is  outlined  below. 

Referrino  to  Fioure  8,  one  may  define  for  the  jth  frame  of  the  ith 

event 


z-.  =  rise  heiaht  from  top  of  bubble  tube 

Dj i  -  diameter  of  center  of  vortex 

nii  =  inclination  of  diameter  to  horizontal 

X°i  i  =  disn^cement  of  center  of  vortex  from  center  of 

bubble  tube 

Dn  initial  bubble  diameter 

'  i 

j  =  1  is  the  first  frame  after  the  bubble  is  broken 

The  following  Physical  parameters  are  averaoed  over  n  events  to 
yield  averaqe  quantities  for  every  jth  frame: 

mean  normalized  rise  height 

mean  normalized  diameter 

mean  rotation 

mean  normalized  translation  of  center  of 
vortex  from  bubble  tube  axis 


19 


Again,  for  a  symmetric  event,  r\.  and  j  should  approach 

In  addition,  standard  deviations,  a  are  computed  for  each 
a  for  some  parameter  P  is  defined  as 


zero  as  a  cneck. 
parameter  where 


1/2 


where 


P. 

J 
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5.  EXPERIMENTAL  RESULTS 


A  discussion  of  laminar  and  turbulent  infinite  fringe  interferometry, 

shadowgraph,  and  particle  tracking  results  are  presented  in  the  folli 
sections. 


lowing 


5.1  HOLOGRAPHIC  INTERFEROMETRY  DATA 

Unlike  the  interferometry  results  reported  in  Section  6.2  of 

Reference  1,  it  should  be  noted  that  all  fringe  distribution  and  helium 

concentration  results  presented  in  this  section  are  based  upon  data 

acquired  using  the  sphere  drop  bubble  bursting  technique  described  in 

Section  2,1  of  this  report  and  the  new  data  analysis  procedure  outlined  in 
Section  4. 1 . 


5-1.1  Laminar  Results 

A  total  of  35  horizontal  view  holographic  interferograns  were 
obtained  of  separate  repeat  events  at  a  tank  pressure  of  one  atmosphere, 

0.19  seconds  after  bubble  bursting.  This  tank  pressure  and  the  initial 
bubble  diameter  of  0.9  inch  correspond  to  a  Reynolds  number  of  1100.  Two 

typical  interferograms  from  this  set  are  presented  for  comparison  purooses 
in  Figures  9  and  10. 

Averaged  radial  fringe  number  profiles  for  left  and  right  sides 
(using  the  method  outlined  in  Section  4)  are  presented  in  Figure  11  and  very 
nearly  coincide,  indicative  of  a  very  symmetric  average  event.  This  average 
event  had  a  rise  of  3.1  initial  bubble  diameters  and  a  diameter  of  2.22 
initial  bubble  diameters. 

Mean  helium  concentrations  corresponding  to  these  data  are  presented 
in  Figure  12.  Since  helium  concentration  is  dependent  on  the  gradients  of 
fringe  number,  wider  discrepancies  between  left  and  right  sides  exist  in 
the  concentration  profiles  than  in  the  fringe  profiles,  though  they  are 
still  reasonably  close.  Note  that  the  peak  concentration  of  about  4  percent 
helium  occurs  at  slightly  less  than  half  of  the  outer  radius. 
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Radial  helium  concentrations  along  cuts  through  the  vortex  wake  were 
calculated  for  the  one  atmosphere  set  of  data  at  1/2,  3/4  and  1  vortex 
diameter  below  the  wake  coordinate  origin  using  the  technioue  described  in 
Section  4.1.4.  The  resulting  mean  concentration  profiles  computed  from 
averaged  fringe  distributions  are  presented  in  Figure  13.  Peak  helium 
concentrations  of  approximately  2.5  percent  exist  near  the  center  of  the 
wake,  become  nearly  a  constant  value  of  about  1.5  percent  through  most  of 
the  wake  and  reduce  to  zero  rapidly  at  the  edge.  It  is  interesting  to 
note  that  peak  helium  concentrations  occurring  at  the  wake  axis  are  about 
half  as  large  as  the  oeak  helium  concentration  within  the  vortex  itself. 

Thus,  even  though  only  a  small  portion  of  helium  forms  the  wake  (the 

wake  volume  is  far  less  than  that  of  the  vortex),  it  remains  quite  concentrated 

5- 1 • 2  Turbulent  Results 

Thirty-three  horizontal  view  holographic  interferograms  were  acquired 
at  a  tank  pressure  of  8  atmospheres,  1.44  seconds  after  bubble  bursting. 

The  Reynolds  number  for  these  events  was  10,400  based  on  this  tank  pressure 
and  ar  initial  bubble  diameter  of  one  inch.  Two  typical  interferograms  from 
this  series  are  pictured  in  Figures  14  and  15. 

Averaged  fringe  numbers  vs.  radius  for  these  data  appear  in  Figure  16 
The  degree  of  agreement  between  the  left  and  right  sides  is  less  than  that 
of  the  one  atmosphere  data,  but  nonetheless  reasonably  close.  This  average 
eight  atmosphere  vortex  had  a  rise  and  diameter  of  9.4  and  4.84  initial 
bubble  diameter,  respectively.  The  aforementioned  limitations  in  the 
infinite  fringe  interferometry  system  prevented  any  complete  analysis  of  a 
high  pressure  vortex  at  an  earlier  stage  in  development. 

Mean  helium  concentrations  as  a  function  of  normalized  radius  were 
calculated  from  these  fringe  distributions  and  are  presented  in  Figure  17 
The  late  stage  in  development  of  the  vortex  is  reflected  by  the  concentration 
profiles  indicating  little  or  no  helium  present  within  the  central  region 
approximately  0.4  times  the  outer  radius  of  the  vortex.  The  size, 
and  the  low  peak  helium  concentration  (compared  to  the  one  atmosphere  case) 
are  indicative  of  the  late  stage  of  development  and  resultant  large  volume 
of  air  entrained  by  the  vortex. 
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5.2  SHADOWGRAPH  SIZE/RISE  STATISTICS 

rate  o/mT  ^  “'V'  movie5  were  obtained  at  a  frame 

rate  of  24  frames/second  at  pressures  of  1,  4,  8,  and  10  atmospheres  to 

etemine  mean  and  standard  deviation  size/rise  characteristics  of  simulated 
e  alls.  They  were  averaged  frame  by  frame  using  the  procedure  outlined 

°"  4-3'  Se,ected  frM,es  *  typical  turbulent  shadowgraph 
mo/le  were  previously  presented  in  Reference  1. 

Non-dimensional  plots  of  rise  vs.  size  are  depicted  on  a  logarithmic 
icale  ,n  Figure  18  for  all  four  Pressures.  Additional  interferometry  data 
ncnnts  and  fireball  size/rise  data  from  an  atmospheric  nuclear  test  are 
else  presented  in  this  figure.  Since  fireball  size  and  rise  are  among 
eMer  lamented  data  from  full  scale  nuclear  events,  an  indication 
o.  the  validity  of  this  laboratory  simulation  may  be  inferred  through  a 
comparison  with  ful,  scale  data.  It  is  apparent  that  the  turbulent  (eight 
and  ten  atmosphere)  laboratory  data  are  in  excellent  ogreement  with  published 
results  of  at  least  one  low  yield,  low  altitude  atmospheric  nuclear  test, 
onversely,  laminar  and  intermediate  Reynolds  number  simulations  reflect 

a  much  faster  rise,  resulting  from  a  less  pronounced  entrainment  of  air 
than  tne  turbulent  events. 

5.3  PARTICLE  TRACKING  RESULTS 

Typical  photographs  of  evolving  particulate  seeded  vortices  are 
presented  in  Figures  19  and  20  for  comparison  purposes  and  illustrate 
dramatically  the  effect  of  turbulence  on  toroid  structure.  Figure  19 
represents  a  laminar  vortex  taken  at  a  pressure  of  one  atmosphere  with  an 
ini  la  bubble  diameter  of  one  inch  corresponding  to  a  Reynolds  number 
of  approximately  1300.  The  important  feature  of  this  photograph  is  the 
distinct  toroid  "jelly-roll"  structure  apparent  throughout  the  entire  rise. 

Figure  20  displays  a  turbulent  event  acquired  at  a  pressure  of  eight 
atmospheres  with  the  remaining  test  Parameters  identical  to  the  one  atmosphere 
hot  graph.  The  Reynolds  number  for  this  particular  event  was  about  10,400 
n  this  case,  turbulent  mixing  has  completely  overwnelmed  all  internal  struc¬ 
ture  evident  in  the  laminar  event.  It  is  apparent  from  these  and  similar 
vortex  structure  photographs  that  a  "jelly-roll”  structure  as  such  is  strictly 
associated  with  a  laminar  event  and  will  not  exist  in  a  full  scale  detonation 
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Figure  12.  Averaged  Radi 


I  oo 


32 


ok 


- 


Figure  16.  Averaged  Radial  Fringe  Distributions  -  Turbulent  Data 
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•  DQ  =  1  INCH 

•  P  =  1  ATMOSPHERE 

•  Re  ^  1300 

•  MICROSPHERE  PART^tS 

•  *5  FLASHES/SEC 


Figure  19.  Multiple  Exposure  Photograph  of  Particulate 
Seeded  Laminar  Fireball 
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•  D  •  1  INCH 

•  MICROSPHERE  PARTICLES 
t  -\-5  FLASHES/SEC 

t  P  =  8  ATMOSPHERE 

•  Re  ^  10,400 


Figure  20.  Multiple  Exposure  Photograph  of  Particulate 
Seeded  Turbulent  Fireball 
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6.  CURRENT  STATUS  AND  FUTURE  PLANS 


Within  the  fringe  spacing  capability  of  the  finite  fringe 
interferometer,  5  to  6  additional  sets  of  i nterferometry  data  will  be 
acouired  and  analyzed  by  the  end  of  this  pfiase  of  the  study.  They  will 
consist  of 

(a)  Minimum  obtainable  Reynolds  number  tests  at  low  rise 

height  and  pressure  to  simulate  the  most  laminar  condition 

(b'  Maximum  obtainable  Reynolds  number  tests  using  N2  bubbles 
in  an  SFg  environment  at  ten  atmospheres  pressure 

(c)  Tests  conducted  at  constant  rise  positions  but  different 
pressures  to  provide  Reynolds  number  dependence  information 

These  tests  are  currently  in  progress  and  results  will  be  published  in  the 
final  technical  report  during  September  1973. 

Direct  turbulence  measurements  using  the  present  diagnostic  instru¬ 
mentation  are  not  possible,  but  an  extension  of  the  study  is  planned  which 
will  allow  simultaneous  velocity  fluctuation  and  density  fluctuation 
measurements  to  be  made  using  a  two-component  laser  Doppler  velocimeter 
and  a  density  probe,  respectively.  These  essential  measurements  will 
be  used  to  guide  and  provide  inputs  for  the  development  of  turbulence 
model  equations  for  fireball  flow  fields. 
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